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f~-«. Abstract 



This paper presents a new numerical approach to the study of non- 
periodicity in signals, which can complement the maximal Lyapunov ex- 
P^ ponent method for determining chaos transitions of a given dynamical 

y;P system. The proposed technique is based on the continuous wavelet trans- 

form and the wavelet multiresolution analysis. A new parameter, the scale 
index, is introduced and interpreted as a measure of the degree of the sig- 
nal's non-periodicity. This methodology is successfully applied to three 
classical dynamical systems: the Bonhoeffer-van der Pol oscillator, the 
logistic map, and the Henon map. 
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1 Introduction 

In the study of chaotic dynamical systems it is quite common to have bifurcation 
diagrams that represent, for each value of one or more parameters, the number 
of periodic orbits of the system. The determination of the parameter values for 
which the system becomes chaotic is a classical problem within the theory of 
dynamical systems [T]. 



Although there is no universally accepted definition of chaos, usually, a 
bounded signal is considered chaotic if (see [2]) 

(a) it shows sensitive dependence on the initial conditions, and 

(bl) it is non-periodic, or 

(b2) it does not converge to a periodic orbit. 

Usually, chaos transitions in bifurcation diagrams are numerically detected 
by means of the Maximal Lyapunov Exponent (MLE). Roughly speaking, Lya- 
punov exponents characterize the rate of separation of initially nearby orbits 
and a system is thus considered chaotic if the MLE is positive. Therefore, the 
MLE technique is one that focuses on the sensitivity to initial conditions, in 
other words, on criterion (a). 

As to criteria (bl) and (b2), Fourier analysis can be used in order to study 
non-periodicity. However chaotic signals may be highly non-stationary, which 
makes wavelets more suitable [S]. Moreover, compactly supported wavelets are 
a useful tool in the analysis of non-periodicity in compactly supported signals, 
as we will see in Corollary [2] 

Wavelet theory is a quite recent area of mathematical research that has been 
applied to a wide range of physical and engineering problems (see, for instance 
[HH] for classical applications to image processing and to time series analysis and 
[H1I71IH] for examples of non-standard applications to DNA sequences, astronomy 
and climatology). In particular, the wavelet decomposition of a signal has been 
proved to be a very useful tool in the study of chaotic systems. Indeed, wavelets 
have been successfully used in the analysis of the chaotic regimes of the Duffing 
oscillator, focusing on the detection of periodicities within chaotic signals, and 
on chaos numerical control (see EH] ) • 

In this paper we present a method for studying non-periodicity that can 
complement the MLE for determining chaos transitions of a given dynamical 
system which can be either discrete or continuous. This method is based on 
the Continuous Wavelet Transform (CWT) and the wavelet Multiresolution 
Analysis (MRA) of a signal. In particular, we compute the ratio of the scalogram 
value at the dominant scale (i.e. the scale where the maximum is reached) to 
the value at the least significant scale (i.e. the scale where the scalogram takes 
its minimum value after that maximum is reached). This quotient determines 
the scale index which is strictly positive when the signal is non-periodic, and 
can be interpreted as a measure of the degree of non-periodicity. 

The paper is organized as follows: Section [2] is devoted to the establishing 
of the background and the main results of the wavelet theory used later in the 
paper. In Section [3] we give an overview of the scalogram and define the scale 
index. Finally, Section |4] illustrates the value of the method by showing the 
successful! detection of the chaos transitions for three classical systems: the 
Bonhocffer-van der Pol (BvP) oscillator, the logistic map and the Henon map. 



2 Wavelet analysis of time series 

In this section we introduce the wavelet tools and results needed for defining 
the scale index. 

2.1 Continuous wavelet transform and scalogram 

Wavelet theory is based on the existence of two special functions 4> and ijj, known 
as the scaling and wavelet functions respectively [TTj . 

A wavelet function (or wavelet, for short) is a function ip ^ L'^ (M) with zero 
average (i.e. LV' = 0)) with \\ip\\ = 1, and centered in the neighborhood of 
i = [TT]. Moreover, we are going to demand that tip (t) E L^ (M) in order to 
ensure that the continuous wavelet transform (l3| is invertible in some way. 

Given a wavelet ip, its dilated and translated dyadic version is given by 

^,.fe(i):=^V'(^^— j, (1) 

where j, fc G Z. It is important to construct wavelets such that the family of 
dyadic wavelets {ipj,k}j,k£Z is an orthonormal basis of L^ (K). These orthonor- 
mal bases are related to the Multiresolution Analysis (MRA) of signals. 

Scaling -0 by a positive quantity s, and translating it by u G K, we define a 
family of time-frequency atoms, ipu,!, as follows: 

ipuAt) ■■= -i'A (^) ' u e M, s > 0. (2) 

Note that there is an abuse of notation in expressions ^ and (IT]) . 

Given f G L^ (M), the continuous wavelet transform (CWT) of / at time u 
and scale s is defined as 

/ + 00 
/(i)v:..(t)dt, (3) 

-oo 

and it provides the frequency component (or details) of / corresponding to the 
scale s and time location t. 

The wavelet transform given in ([3]) , provides a time-frequency decomposition 
of / in the so called time-frequency plane (see Figure fl]). 

The scalogram of /, S, is defined as follows: 

S is) := \\Wf (u, s)\\^( f " \Wf {u, s) \^du 



S (s) is the energy of the continuous wavelet transform of / at scale s. Obviously, 
S{s) > for all scale s, and if S{s) > we will say that the signal / has details 
at scale s. Thus, the scalogram is a useful tool for studying a signal, since it 
allows the detection of its most representative scales (or frequencies), that is, 
the scales that mostly contribute to the total energy of the signal. 
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Figure 1: Time-frequency plane decomposition corresponding to the BvP solu- 
tion with A = 0.76 (left) and A — 1.16 (right) using Daubechies (eight-wavelet 
and four-wavelet respectively) wavelet functions (see Section 4.1 for the defi- 



nition of the BvP system). Each point in this 2D representation corresponds 
to the modulus of the wavelet coefficients of the CWT. Note that the wavelet 
coefficients of the CWT with A = 1.16 vanish at scale 12 (i.e. twice its period) 
at any time, as we will prove in Theorem IT] 



2.2 Analysis of compactly supported discrete signals 

In practice, to make a signal / suitable for a numerical study, we have to 
(i) consider that it is defined over a finite time interval / — [a, 6], and 
(ii) sample it to get a discrete set of data. 

Regarding the first point, boundary problems arise if the support of ^u,s 
overlaps t — a or t = b. There are several methods for avoiding these problems, 
like using periodic wavelets, folded wavelets or boundary wavelets (see [H]); how- 
ever, these methods either produce large amplitude coefficients at the boundary 
or complicate the calculations. So, if the wavelet function ip is compactly sup- 
ported and the interval I is big enough, the simplest solution is to study only 
those wavelet coefficients that are not affected by boundary effects. 

Taking into account the considerations mentioned above, the inner scalogram 
of / at a scale s is defined by 



S" 



(s):=||Ty/(s,u) !!,,(,) = 




\Wfis,u)\Mu 



where J{s) — [c{s),d{s)] C / is the maximal subinterval in I for which the 
support of ipu^s is included in / for all u € J{s). Obviously, the length of / must 
be big enough for J{s) not to be empty or too small, i.e. b — a^ si, where I is 
the length of the support of ip. 



Since the length of J{s) depends on the scale s, the values of the inner 
scalogram at different scales cannot be compared. To avoid this problem, we 
can normalize the inner scalogram: 



cinncr / ^ 

S (s) - 



-inner 



(d(s)-c(s))i 



With respect to the sampling of the signal, any discrete signal can be an- 
alyzed in a continuous way using a piecewise constant interpolation. In this 
way, the CWT provides a scalogram with a better resolution than the Discrete 
Wavelet Transform (DWT), that considers dyadic levels instead of continuous 
scales (see [TT]). 

3 The scale index 

In this section we introduce a new parameter, the scale index, that will give us 
information about the degree of non-periodicity of a signal. To this end we will 
first state some results for the wavelet analysis of periodic functions (for further 
reading please refer to [TT] and references therein). 

If / : K — >■ C is a T-periodic function in L^ ([0, T]), and -0 is a compactly sup- 
ported wavelet, then Wf {u, s) is well-defined for u S R and s G K"*", although 
/ is not inL'^iR). 

The next theorem gives us a criterion for distinguishing between periodic 
and non-periodic signals. It ensures that if a signal / has details at every scale 
(i.e. the scalogram of / does not vanish at any scale), then it is non-periodic. 

Theorem 1 Let / : M — > C be a T-periodic function in L^ ([0, T]), and let Tp be 
a compactly supported wavelet. Then W f {u, 2T) = for all u € M. 

For a detailed proof see Appendix|A] From this result we obtain the following 
corollary. 

Corollary 2 Let f : I — [a,b] ^ C a T-periodic function in L^ ([a, a -\- T]). If 
ip is a compactly supported wavelet, then the (normalized) inner scalogram of f 
at scale 2T is zero. 

These results constitute a valuable tool for detecting periodic and non- 
periodic signals, because a signal with details at every scale must be non-periodic 
(see Figure l2|. Note that in order to detect numerically wether a signal tends to 
be periodic, we have to analyze its scalogram throughout a relatively wide time 
range. 

Moreover, since the scalogram of a T-periodic signal vanishes at all 2kT 
scales (for all k E N), it is sufficient to analyze only scales greater than a 
fundamental scale sq. Thus, a signal which has details at an arbitrarily large 
scale is non-periodic. 

In practice, we shall only study the scalogram on a finite interval [so,si]. 
The most representative scale of a signal / will be the scale Smax for which the 



A=0.60(l-periodic) 



A=0.66 (2-periodic) 




Figure 2: Normalized inner scalograms for certain solutions of the BvP system 



(Section |4JJ), from t = 20 to t = 400 {At = 0.05), for different values of A, the 
scale parameter s running from Sq = 0.05 to si = 12.8, with As = 0.05, and 
using Daubechies eight-wavelet function. It is observed how the scalogram of 
T-periodic signals vanishes at s = 2T. 



scalogram reaches its maximum value. If the scalogram S{s) never becomes 
too small compared to 5(smax) for s > Smax, then the signal is "numerically 
non-periodic" in [so,si]. 

Taking into account these considerations, we will define the scale index of / 
in the scale interval [sq, Si] as the quotient 



^scalc ■ 



S[Saun) 
^ l.'^niax j 



where Smax is the smallest scale such that S{s) < 5(s„iax) for all s G [so,si], 
and Sinin the smallest scale such that 5(sniin) < '5(s) for all s G [smax, si]- Note 
that for compactly supported signals only the normalized inner scalogram will 
be considered. 

From its definition, the scale index igcaic is such that < iscaic < 1 and it 
can be interpreted as a measure of the degree of non-periodicity of the signal: 
the scale index will be zero (or numerically close to zero) for periodic signals 
and close to one for highly non-periodic signals. 

The selection of the scale interval [sq, Si] is an important issue in the scalo- 
gram analysis. Since the non-periodic character of a signal is given by its be- 
havior at large scales, there is no need for sq to be very small. In general, we 
can choose sq such that Smax = sq + e where e is positive and close to zero. 
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Figure 3: Left: normalized inner scalogram of the almost-periodic function 
sin(i) + sin (t/V^) + sin (t/V5), from i = to i = 800 (Ai = 0.1), using the 
Daubechies four-wavelet function. Right: the scale index tends to zero as we 
increase the scale si. 



On the other hand, si should be large enough for detecting periodicities. 
For example, if we have an almost-periodic function / defined on M, then for 
any given e > there exists an almost-period T (e) such that 

\f{t + T)-f{t)\<e 

for all t e M (see [12]). Hence, if we choose e numerically close to zero, there 
is a large enough value of T for which the function is "numerically T-periodic" , 



and the scale index will be close to zero if si is greater than 2T. So, for an 
almost-periodic function, the scale index tends to zero as we increase si (see 
Figure [3]). But as si increases, so does the computational cost. In fact, the 
larger si is, the wider the time span should be where the signal is analyzed, in 
order to maintain the accuracy of the normalized inner scalogram. 

Scales Smin and Smax determine the pattern that the scalogram follows (see 
Figure p| . For example, in non-periodic signals Smin can be regarded as the 
"least non-periodic scale". Moreover, if s^nin — si, then the scalogram decreases 
at large scales and si should be increased in order to distinguish between a 
non-periodic signal and a periodic signal with a very large period. 




Figure 4: Values of Smax (grey) and Smin (black) for orbits of the BvP oscillator. 



4 Examples 

In this section we illustrate how the scale index igcaio is used in order to detect 
and study non-periodic orbits of three classical dynamical systems: the logistic 
map, the Henon map and the forced Bonhoeffer-van der Pol oscillator. 

The reason for choosing these dynamical systems as examples for testing 
the validity of the scale index is mainly that they are three well known chaotic 
dynamical systems, arising from different research areas, and are mathematically 
very different. These systems present typical bifurcation diagrams with chaotic 
and non-chaotic regions. In order to show the effectiveness of the index iscaio, 
we compare the bifurcation diagram, the MLE, and iscaio- It will be shown 
that there is a correspondence between the chaotic regions of the bifurcation 
diagram, the regions where the MLE is positive, and the regions where iscaio is 
positive. 

Figures [s] and [6] depict the comparison between the three methods mentioned 
above. The signals were studied from Iq — 20 to identify not only periodic 
signals, but also signals that converge to a periodic one. For the computation of 
the MLE, 1500 iterations were used. Integer scales between sq = 1 and si = 64 
were considered in the computation of iscaio- A scale was considered to have no 
details if the scalogram at that scale takes a value below e = lO"'*. 



4.1 Continuous Dynamical System: The Bonhoeffer-van 
der Pol Oscillator 

The BonhoefFer-van der Pol oscillator (BvP) is the non-autonomous planar sys- 
tem 

x' = x-^-y + I{t) 
y' = c{x + a — by) 

being a, b, c real parameters, and / (t) an external force. We shall consider 
a periodic force I{t) ~ A cos (27rt) and the specific values for the parameters 
a = 0.7, b = 0.8, c — 0.1. These values were considered in [T5] because of their 
physical and biological importance (see [Ml). 

The classical analysis of the BvP system is focused on its Poincare map, 
defined by the flow of the system on t = 1 (see [T]). Plotting the first coordi- 
nate of the periodic fixed points of the Poincare map versus the parameter A 
(amplitude of the external force), a bifurcation diagram is obtained |15| . 

Such diagrams present chaotic and non-chaotic zones. From a geometric 
point of view, chaos transitions are related to homoclinic orbits (creation or 
destruction of Smale horseshoes) between the invariant manifolds of a saddle 
fixed point of the Poincare map (see [T]). Such a relationship is thoroughly 
described in a recent work [IB] . 

Figure [5] includes the analysis of the BvP system. The parameter range has 
been split in two regions, 0.7 < A < 0.8 and 1 < ^ < 1.3, which are the regions 
were chaotic orbits are found. Note the high level of agreement between the 
MLE and the igcaie index: the values of A for which the MLE is negative are 
also the values for which igcaie ~ 0. 

Also remarkable is the coincidence between a relative maximum in the Zgcaie 
and the well known "sudden" expansion of the size of the attractor at ^4 » 
0.748 (see [16], [13]). Moreover, another relative maximum fits with a "sudden" 
contraction at A^i 1.27. 

4.2 Discrete Dynamical Systems: The Logistic and the 
Henon maps 

The logistic map is the discrete dynamical system given by the difference equa- 
tion 

xt+i = Axt{l - xt). 

This well known dynamical system, arising from population dynamics theory, 
is a classical example of a simple polynomial map whose orbits exhibit chaotic 
behavior for some values of the real parameter A. 

On the other hand, the Henon map is the two-dimensional dynamical system 
given by the quadratic map 

xt+i = 1~ Axl + yt 
yt+i = bxt 
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Figure 5: Comparison between the bifurcation diagram, MLE and Zgcaio (from 
top to bottom) for the BvP osciUator. 
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with A and b real parameters. For the values A = 1.4 and b — 0.3, giving what 
is often called the canonical Henon map, a strange attractor is present. Fixing 
the value b = 0.3 and varying the parameter A, the map may be chaotic or not. 
Through the representation of the x-coordinate of the orbits versus the value of 
the parameter A, a typical bifurcation diagram is obtained. 

Figure [6] shows the comparison between the bifurcation diagram, the MLE 
and iscaie for both the logistic map (left) and the Henon map (right). As in the 
case of the BvP oscillator, the agreement between the MLE and igcaio is clear. 

It is also noticeable that the maximum values of the iscaio in these two 
dynamical systems are reached when the main branches of their bifurcation 
diagrams overlap. 
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Figure 6: Comparison between the bifurcation diagram, MLE and iscaic (from 
top to bottom) for the logistic map (left) and the Henon map (right). 
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5 Conclusions 

Wavelet analysis has proved to be a valuable tool in the study of chaotic systems. 
In particular, scalogram analysis and the introduced scale index are a good 
complement to the MLE, since the scale index gives a measure of the degree 
of non-periodicity of the signal, while the MLE gives information about the 
sensitivity to initial conditions. Thus, the combination of the two methods 
gives us a comprehensive description of a chaotic signal. 

Since the MLE and the scale index focus on different characteristics of the 
signals, the latter contributes to detecting effects that are not detected by the 
MLE, for example, the sudden expansion of the size of the attractor in the BvP 
system, and the overlapping of the main branches of the bifurcation diagram 
in the logistic and Henon maps. Moreover, there are regions in the logistic and 
Henon graphs where the MLE is increasing but the iscaic is decreasing. This 
means that while the sensitivity to the initial conditions is greater, on the other 
hand the signal is less non-periodic. 

Additionally, the study of the scale index does not require an analytical 
expression for the signal. In those cases where an analytical description of the 
dynamical system involved is not available (e.g. experimental signals), although 
there are methods for estimating the MLE (see [17]), the scale index might be 
a useful alternative. 

These techniques can also be applied to any other discipline where the anal- 
ysis of time series is required, such as Earth sciences, econometry, biomedicine 
and any other one where non-linear behavior is expected to occur. For instance, 
this opens a branch to the study of colour noise presented in this kind of time se- 
ries; it also could be used to determine properties of seismic and volcanic events 
f|18|): to detect chaos in non- linear economic time series (PJj) as a complement 
to other methods already used, like those based on Lyapunov exponents and the 
power spectral density; and in biomedicine, where wavelets have already been 
used to analyze non-stationary cardiac signals f|20p. 

A Proof of Theorem [1] 

Given g : M — > C with compact support, its periodization over [0, 1] is defined 
as 

gP-(i):^^g(t-ffc), teR. 

feez 

It is clear that, if ^^(t) '.— g{t — u) is a translated version of g, then 

{gur'{t)^g''''it-u). (4) 

Moreover, if /i e L^ ([0, 1]), then 

(/jP-,5) = (;i,5P-)[o,i]. (5) 

Lemma 3 The periodization z/'ST '-^ ■^c™ almost everywhere for u = and 
s = 2. 
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Proof. In this proof we are going to use the dyadic notation for the dilated and 
translated wavelets given by expression (HI). So, we have to prove that -0^'^'^ = 
a.e. for j = 1 and fc = 0. 

Let /i £ i^ (fO, ll). Since \ ib^-f. \ Ujlj is an orthogonal basis 

I •''''J-oo<i<0.0</c<2-J 

of L^ ([0, 1]) and the periodization ip^^^ is orthogonal to the family of functions 

(see [m Thm. 7.16, Lem. 7.2]), we have 



{€r}. 



oo<j<0, 0</c<2-J 

(/i,^P^J-)[o,i] =M(l,V'ro)[o,i] = Af(l,^i.o), 



where M — (ft,, l)ro.ii. Taking into account that tp has zero average, we have 
that (1,^1 o) = 0, and so 

(^, Co) [0,11=0. 

In general, it can be proved that -0''" = a.e. for all j > 1 and for all fc € Z. 



Lemma 4 The periodization "4)^2 — ^ o,-^- /<"" aH u G M. 

Proof. Since V'„^2(^) ~ "00. 2(^ — "") we have that i^^l{t) = V-'o'2 (* — m) = a.e. 
taking into account Q and Lemma Is] 

■ 

Theorem fll Let f -.R^ C be a T-periodic function in L^ ([0,T]), and let 
ip be a compactly supported wavelet. Then W f (u, 2T) = for all m G K. 
Proof. Without loss of generality wc may assume that T = L Then, it will be 
proved that Wf {u, 2) = for ah u g E: 

Let h := /|[o,i]- Then / = h^'^'^ and we get 

Wf{u,2) = (/,0„,2> = (/iP°^^u,2> = {h,ru%)[o,i] = 0, 
taking into account ([5]) and Lemma HI 
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